Temporal and diffraction effects in entanglement creation in an optical cavity 
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A practical scheme for entanglement creation between distant atoms located inside a single-mode 
optical cavity is discussed. We show that the degree of entanglement and the time it takes for the 
entanglement to reach its optimum value is a sensitive function the initial conditions and the position 
of the atoms inside the cavity mode. It is found that the entangled properties of the two atoms 
can readily be extracted from dynamics of a simple two-level system. Effectively, we engineer two 
coupled qubits whose the dynamics are analogous to that of a driven single two- level system. It is 
found that spatial variations of the coupling constants actually help to create transient entanglement 
which may appear on the time scale much longer than that predicted for the case of equal coupling 
constants. When the atoms are initially prepared in an entangled state, they may remain entangled 
for all times. We also find that the entanglement exhibits an interesting phenomenon of diffraction 
when the the atoms are located between the nodes and antinodes of the cavity mode. The diffraction 
pattern of the entanglement varies with time and we explain this effect in terms of the quantum 
property of complementarity, which is manifested as a tradeoff between the knowledge of energy of 
the exchanged photon versus the evolution time of the system. 

PACS numbers: 32.80.Qk, 42.50.Fx, 42.50.Dv 
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I. INTRODUCTION 



The study of practical schemes for creation of quantum 
entanglement between atoms (or ions) is the most active 
area in quantum optics and quantum information sci- 
ence [l| . Different schemes have been proposed including 
atom trapped inside a single mode cavity v & b i Ul, or in- 
side two separate cavities d, H, 0, H, H, fltxTllL . One 
of the most popular scheme involves two two-level atoms 
located (trapped) within a single-mode cavity field, ft 
has been demonstrated that the entanglement could in 
principle be created through a continuous observation of 
the cavity field 11311 or through dispersive atom-cavity 
field interactions [l4, ITU [HI], thereby creating a strong 
" action at a distance" . The approach used is straightfor- 
ward: Provided no photon is leaking through the cavity 
mirrors or no photon is exchanged between the atoms 
and the cavity field, a pure entangled state between the 
two atoms results. However, all of these procedures for 
generating entangled atoms have suffered from a com- 
mon handicap: their choice of equal coupling strengths 
of the atoms to the cavity mode. The difficulty is that the 
entanglement depends on the the coupling constant be- 
tween the atoms and cavity mode which depends, in turn, 
on the location of the atoms in the cavity mode. In a 
standing wave cavity, one can achieve the equal coupling 
constants by locating the atoms precisely at the antin- 
odes of the cavity field, or by sending slowed (cooled) 
atoms through a cavity antinode in the direction perpen- 
dicular to the cavity axis. This is a relatively easy task 
at microwave frequencies and, in fact, detection of entan- 



gled atoms have already been performed on a beam of 
Rydberg atoms traversing a superconducting microwave 
cavity [17] . However, at optical frequencies this task may 
well be hard to achieve. In practice, the dipole coupling 
constants vary with the location of the atoms in the cav- 
ity mode. For example, in a standing-wave structure of 
the cavity mode, the coupling constant varies with the 
position of the atom as [l8| 



_ 2 / 2 

g(r) = g(r, z) = goe r /w ° cos (kz) 



(I) 
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where z determines location of the atom along the cav- 
ity axis, k = 2tt/\ is the wave number of the field, 
r = (x 2 + j/ 2 ) 1 / 2 is the distance of the atom from the 
cavity axis, and wq is the mode waist. The coupling 
constant reaches maximum value go when the atom is lo- 
cated on the cavity axis at an antinode of the standing 
wave. In practice one would like the varying coupling 
constant g(r) to coincide with go. However, because of 
the small wavelength, locating the atoms precisely at the 
antinodes of the cavity mode and thereby eliminating the 
variation of g(r) with position is very difficult in an opti- 
cal cavity. This may change the physics completely and 
thus suggests that the problem of creating entanglement 
between distant atoms in an optical cavity is a significant 
experimental challenge. 

In this paper, we are concerned with the optical frequency 
regime and investigate the mechanism involved in cre- 
ation of entanglement between distant atoms coupled to 
a single-mode cavity field. Our treatment closely follows 
the approach that was used Refs [3, EE EH , but with 
one essential difference. We include a possible variation 
of the coupling constant g(r) with the location of the 
atoms in a standing-wave cavity mode. We are particu- 
larly interested in the consequences of this variation on 
entanglement creation between the atoms, since this will 
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be very pertinent to any practical experimental arrange- 
ments as the distances involved are very small. First of 
all, we derive the general master equation for the reduced 
density operator p s of two two-level atoms coupled to a 
single mode cavity field. Our approach holds for atoms at 
rest or slowly moving through a single-mode cavity. For 
sufficiently cold or slow atoms, radiative equilibrium is 
reached with an essentially fixed coupling constant, g(f), 
at every point inside the cavity, permitting studies of en- 
tangled properties of the system without performing the 
average over random locations and atomic dipole orienta- 
tions. We assume the atoms are far enough apart that the 
direct dipole-dipolc coupling or other direct interactions 
between the atoms can be neglected. This also allows 
a selective preparation of the atoms such that a given 
set of initial conditions for the atomic states is achieved. 
We solve the master equation for two atoms coupled to a 
cavity mode in the limit of a large detuning of the cavity 
mode from the atomic resonance frequency. This enables 
to eliminate the cavity mode and obtain a new master 
equation in which we will recognize some terms equiva- 
lent to the dipole-dipole interaction between the atoms 
and to the Stark shift of the atomic resonance frequen- 
cies. The analogy of this system with that of a single 
two-level atom driven by a detuned coherent field is ex- 
ploited and discussed. This analogy provides a simple 
description of the process of entanglement creation and 
leads to a useful pictorial representation of the system 
in terms of the Bloch vector model. To quantify the 
degree of entanglement, we use the concurrence that is 
the widely accepted measure of two-atom entanglement. 
Simple analytical expressions are obtained for the concur- 
rence that are valid for arbitrary initial conditions and 
arbitrary positions of the atoms inside the cavity mode. 
We obtain the interesting result that spatial variations of 
the coupling constants actually help to create transient 
entanglement which may appear on the time scale much 
longer than that predicted for the case of equal coupling 
constants. We explain this effect in terms of the degree of 
localization of the energy induced in the field by the in- 
teracting atoms. Moreover, we find that for an imperfect 
location of the atoms inside the cavity mode the entangle- 
ment exhibits an interesting time-dependent diffraction 
phenomenon. 



II. MASTER EQUATION 

The system we consider consists of two identical two- level 
atoms (qubits) with upper levels | } , (i = 1,2), lower 
levels |<7i), and separated by energy ftw®. The atoms are 
coupled to a standing-wave cavity mode with the position 
dependent coupling constants g(r{), and damped at the 
rate 7 by spontaneous emission to modes other than the 
privileged cavity mode. The cavity mode is damped with 
the rate k and its frequency lo c is significantly detuned 
from the atomic transition frequency u>o, so there is no 
direct exchange of photons between the atoms and the 



cavity mode. The behavior of the system is described by 
the density operator p, which in the interaction picture 
satisfies the master equation 

| = 4[/M + ~7A^+~^cP, (2) 

where 

2 

H = hJ2[9(r ] )aS+e- tAt +R.c] (3) 

describes the interaction between the cavity field and the 
atoms, 

C c p = (2apa) — a) ap — pa) a) , (4) 

and 

2 

C aP = J2 (2SjpS+ - S+Sjp - pS+Sj) . (5) 

are operators representing the damping of the atoms by 
spontaneous emission and of the cavity field by the cav- 
ity decay, respectively. The operators and S~ are 
the raising and lowering operators of the jth atom, Sj 

describes its energy, a and a) are the cavity-mode an- 
nihilation and creation operators, A = lu c — u>q is the 
detuning of the cavity-mode frequency from the atomic 
transition frequency, and r) is the position coordinate of 
the jth atom within the cavity mode. 
The atoms located at different positions may experience 
different coupling constants, that is g{f\) ^ 3(^*2)- Be- 
cause it is precisely the effect of unequal coupling con- 
stants that interest us most here, we choose the reference 
frame such that 

g{ri)=go, and g(r 2 ) = go cos (kr 12 ) , (6) 

where r\i — z 2 — z\ is the distance between the atoms. 
This choice of the reference frame corresponds to a sit- 
uation where atom 1 is kept exactly at an antinode of 
the standing wave and the atom 2 is moved through suc- 
cessive nodes and antinodes of the standing wave. This 
choice, of course, involves no loss of generality. 
We also assume that the atoms are stationary during the 
interaction with the cavity mode, i.e. the distance be- 
tween the atoms is independent of time (the Raman-Nath 
approximation). This is a good approximation for many 
experiments on cooling of trapped atoms, where the stor- 
age time of the trapped atoms is long, so that they are 
essentially motionless and lie at known and controllable 
distances from one another fl9| . 

In order to study the dynamics of the system, we intro- 
duce density-matrix elements with respect to the cavity 
field mode, denoting (n\p\m) by p nm , and find from the 
master equation ^ that the populations of the two low- 
est energy levels and coherence between them satisfy the 
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following equations of motion 



where 



Poo 



2 _ 1 



3=1 



poi = iApoi - iy^Qj (Sfpn - poo-S/) 

3=1 

1 r - 1 - 

+ An - 2«Poi, 

2 

Pio = -«Apio + i^^gj (puSj - Sj poo) 

3=1 

1 - - 1 - 

+ 2^0^10 - 2 K Pio> 

2 

pn = ^E*^^ 1- ^ 10 ^) 

3=1 

+\^C a pn - npn, 



(7) 



where gj = <?(Tj) and p nm are slowly varying parts of the 
coherences, poi = poi exp(iAt) and pio = Pio exp(— iAt). 
We have considered only the two lowest energy levels as 
in the limit of a large detuning A only the ground state 
(n = 0) and the one-photon state (n = 1) of the cavity 
mode are populated. 

Now, we explicitly apply the adiabatic approximation 
that for a large detuning, the coherences poi and pio 
vary slowly in time, so we can assume that poi ~ and 
pio ~ 0. In this case, we find from ([7]) that in the limit 
of A > gj > 7, k 



Poi 



Pio 



I 2 



3=1 
2 



A E^ (PuSj - Sj Poo) 



(8) 



3=1 



Knowledge of the coherences poi and pio allows us to de- 
rive the master equation for the reduced density operator 
of the atoms. It is done in the following way. First we 
substitute ([5]) into 0, and after neglecting the popula- 
tion pn, as the cavity mode will never be populated, we 
find 



Poo 



% 1 
= x 9i9j [S+S~,p 00 ] + -jC a poo- (9) 



i,j=l 



Since poo = Ttf(p) = p s is the reduced density opera- 
tor of the atoms, we obtain the master equation for the 
density operator of the atoms 



t=l i^3=l 



+ lp£aPsi 



(10) 



9i_ 

A' 



and f2 7 - 



A ' 



(11) 



The first two terms in the master equation (jlOp depend 
on the position coordinate of the atoms and give rise to 
frequency shifts of the atomic levels and the coupling be- 
tween the atoms, respectively. The third term represents 
the damping of the atoms through the interaction with 
the environment. The parameter Si represents the shift 
in energy separation of the levels of the ith atom due to 
the dispersive interaction with the cavity mode. It is an 
analog of a dynamic Stark shift. One can easily see from 
the structure of the first term in the master equation that 



j=i i=i 

2 

= iJ26 i [S?,p s ], 



(12) 



which clearly shows that this term is an analog of the 
energy shift term. Thus in the interaction picture used 
here, it represents a shift of the atomic energy levels. We 
note that due to nonequivalent positions of the atoms, the 
shift of the energy levels is different for different atoms. 
The multi-atom term fi^ represents the shift in energy 
separation of the levels of atom i due to its interaction 
with the atom j through the cavity mode. If the atoms 
are located at antinodes of the standing wave, the term 
f2y is maximal, whereas fiy = if at least one of the 
atoms is located at a node of the standing wave. From 
the structure of the second term in Eq. (fTU)) one can rec- 
ognize that fiij is an analog of the familiar dipole-dipole 
interaction between the atoms 0, [2l[ • This shows that 
the interaction of the atoms with a detuned cavity fields 
produces a structure in the master equation analogous to 
the dipole-dipole interaction between the atoms. 
The above procedure shows that the adiabatic elimina- 
tion of the cavity mode creates a shift of the atomic tran- 
sition frequencies and an effective interaction between 
two distant atoms. Thus, the dynamics of the system 
composed of two identical atoms in nonequivalent po- 
sitions in the cavity mode is equivalent to those of two 
non-identical atoms of different transition frequencies. In 
other words, the procedure is an example of how one can 
"engineer" the dipole-dipole interaction between distant 
atoms. It is now easy to understand why two indepen- 
dent atoms coupled to a strongly detuned cavity mode 
can exhibit entanglement. Simply, the reduced system 
is equivalent to two atoms coupled through the induced 
dipole-dipole interaction. Exactly, this process gives rise 
to the entanglement. 

We point out in passing that despite of the presence of 
the coherent dipole-dipole interaction term, the master 
equation (|10|) is not fully equivalent to the master equa- 
tion of two collective interacting atoms [23| . This is be- 
cause there is no contribution from the cross-damping 
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terms involving dipole operators of two different atoms. 
In other words, the interaction with a strongly detuned 
cavity mode does not create the collective damping of 
the atoms. As a result, the atoms interact independently 
with the environment, so that the system does not evolve 
to a dark state characteristic of the completely collective 
system (2^, . The cross-damping terms would ap- 

pear in the master equation if one assumes the near res- 
onant interaction, A w 0, and the "bad-cavity" limit of 
«»9j»A|i|. 

III. EQUIVALENT TWO-LEVEL DYNAMICS 

The question we are interested in concerns the conse- 
quences of the spatial variation of the coupling constant 
g on the entanglement creation between two atoms lo- 
cated in a strongly detuned single-mode cavity field. To 
answer this question we consider the evolution of the di- 
agonal density matrix elements which correspond to the 
occupation probabilities of the energy levels of the two- 
atom system. Using the master equation (|10j) . we find 
the following equations of motion 

Pu = 7 - 7 On + Pu) , (13) 

p 4 4 = -27P44, (14) 
f>22 = -7P22 + 7P44 + «^12 (P23 - P32) , (15) 
P33 = -7P33 + 7P44 - «^12 (p23 - P32) , (16) 

while the off-diagonal density matrix elements that are 
coupled to the diagonal elements obey the equations 

f>23 = - (7 _ ^12) P23 + «Ol2 (p22 - P33) , (17) 
P32 = ~ (7 + iSu) P32 - «^12 (P22 - P33) , (18) 

where we have used the standard direct-product basis 
given by 

|1) = \ 9l )\g 2 ), |2> = | 5l )|e 2 ) ! 

|3) = |ei>|sa>, |4) = |ei)|e 3 ). (19) 

Here, 5±2 = <5i — 62 is a difference between the single-atom 
Stark shifts. This parameter is of central importance here 
as it determines the relative variation of atomic transition 
frequencies with position of the atoms inside the cavity 
mode. In the special case of g\ = <?2, the parameter 
<5i2 = 0, but this can happen only when the atoms are in 
equivalent positions inside the mode. 
It is easy to see that the set of the coupled equations (|13[) - 
(|18[) can be split into two independent sets of equations 
of motion, which can be written in the form 

Pu = 7-70n+P44), 

P44 = -27/544, 

P++ = -7P++ + 27P44, (20) 

and 

U = — 7U + (5i2W, 

v = — 7W — S12U — 2O12W, 
w = -jw + 2Q 12 v, (21) 



where p ++ = p 22 + P33, u = p 23 + p 32 , v = i{p 23 - P32), 
and w = p 22 — p 33 . Note that the equations of motion 
(f2"Tj) are the exact equivalent of the optical Bloch equa- 
tions of a two-level system driven by a detuned coherent 
field, where the dipole-dipole interaction Qi2 couples to 
the levels like a coherent field, and S12 appears as a de- 
tuning of the field from the driven transition. The upper 
and lower energy levels |2) and |3) thus show a dynam- 
ics that is analogous to that of a driven two-level sys- 
tem. It should be pointed out that the analogy is not 
absolute because, in contrast to the case of a single two- 
level atom driven by a detuned laser field, fli 2 and ^12 
in Eq. (|21[) depend on the same parameters 31 and 172, 
i.e. on the coupling constants of the atoms to the cavity 
field. Consequently, the parameters fl± 2 and S12 are not 
independent. 

The dynamics of the effective two-level system can be 
easily visualized in the Bloch vector model [25|. In this 
model, the system and the driving field are represented 
by vectors in a three-dimensional space, and the time evo- 
lution is simply visualized as a precession of the system- 
state vector about the driving field. In terms of the Bloch 
vector, Eq. (|2ip can be written as 

^ = -75 + d B x 5, (22) 
at 

where (Ib = {— 2f2i2 5 0) ^12) is the pseudofield vector and 
B = (u, v, w) is the Bloch vector. The quantities u and v 
are, respectively, the real and imaginary parts of the co- 
herence between the levels |2) and |3), and w is the popu- 
lation inversion. In the present problem, the coherence is 
induced by the dipole-dipole interaction Q12 which, as we 
have already pointed out, plays a role similar to the Rabi 
frequency of the coherent interaction between the levels, 
i.e. represents a rate at which one quantum of excitation 
is exchanged between the atoms. It should be noted that 
there is no electric-dipole moment between the levels |2) 
and |3), so there are no radiative transitions between the 
levels of the two-level system. The damping rate 7 that 
appears in Eq. (|2"Tj) represents spontaneous decay out of 
the system to the auxiliary level |1). 
Since the u and v components of the coherence are re- 
lated to the interaction between the atoms, their prop- 
erties should be reflected in the entanglement between 
the atoms. In order to determine the amount of entan- 
glement between the atoms and the entanglement dy- 
namics, we use concurrence that is the widely accepted 
measure of entanglement. The concurrence introduced 
by Wootters [26[ is defined as 

C = max(0, a/aT- \/% - ~ \Ai) , (23) 
where {Aj} are the the eigenvalues of the matrix 

R = PsPs, (24) 

with p s given by 

P S = <Ty <£> CTy P* (Ty ® (Jy, (25) 
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and a y is the Pauli matrix. The range of concurrence 
is from to 1. For unentangled atoms C — whereas 
C = 1 for the maximally entangled atoms. In terms of 
the density matrix elements, the concurrence is given by 

C{t) = 2max{o, \p 23 (t)\ - vW(*R<)} , ( 26 ) 

which shows that the basic dynamical mechanism for 
entanglement creation in this system is the coherence 
P23(t)- That is, the origin of the entanglement in the 
system can be traced back to the time evolution of the 
coherence P2z{t)- 

Utilizing the relation P2z{t) = ( u (t) — the time 

evolution of the coherence p 2 3(t) can be readily found 
from the solution of the Bloch equations (HJ). The gen- 
eral solution for u(t) and v(t), valid for arbitrary initial 
conditions, is given by 

u{t) = ult)^ 1 
e~ 7 * 

= — — [2Qi2-4 + 8\2 (voasinat + 2? cos erf)] , 
or 

e -7 ' 

v(t) = v(t)e^ 7t = [v a cos at + B sin at] , (27) 

a 

where wq = w(0), uq = u(0) and vq = v(0) determine the 
initial population inversion and coherences in the system 

A = 2f2i 2 uo - S12W0, B = 5i 2 u + 2SI12W0, (28) 

and a = ^Afl 2 2 + S 2 2 is the detuned Rabi frequency. 
It is also worthwhile to find the time evolution of the 
population inversion 

w(t) = w(t)e~^ = e —{-5 12 A 

OT 

+2fii2 (v a sin at + B cos at)}. (29) 

which will allow us to use the Bloch vector picture to 
gain insight into the problem of entanglement creation 
in the system. For a single-quantum excitation, pn(t) = 
P44(t) = 0, and then the Bloch vector component w(t) 
determines the population distribution among the levels 
1 2) and |3), and through the relation 

u 2 (t)+v 2 (t)+w 2 (t) =e" 27 *, (30) 

the component, in turn, determines the entanglement 

C{t) = max {0, e-" v / l-w 2 (t)} . (31) 

This implies that the concurrence can be completely de- 
termined by observing changes in the populations of the 
system's energy levels. When the population is in the 
level 1 2) or in the level 1 3), w{t) = ±1, and then C(t) = 0, 
whereas C(i) achieves its optimum value C(t) — 1 when 
w(t) = 0. Therefore, we can interpret the entangle- 
ment as a consequence of a distribution of the population 
among the energy levels. 



In the complete Bloch vector picture, the Bloch vec- 
tor makes a constant angle 9 = tan _1 (— 2f2 12 /<5i2) with 
fls , it rotates around it, tracing out a circle on the 
Bloch sphere. When the Bloch vector does not inter- 
sect the "north pole" B = (0,0,1) or the "south pole" 
B = (0,0,-1), the inversion w(t) ^ ±1 and then an 
entanglement between the atoms occurs. 
It is clear from Eqs. (|29p and (|3"Tj) that the temporal evo- 
lution of the concurrence can be quite different depending 
on the initial conditions. However, the most interesting 
aspects of the solutions (|29|) and (|3Tj) relate to the mod- 
ifications of the time evolution of the system produced 
by the detuning ^12, because these features are not en- 
countered at all under the equal coupling constants and 
appear never to have been investigated before. This mo- 
tivates our study of the effect of a spatial location of the 
atoms on the atom-atom entanglement. 

IV. ENTANGLEMENT DYNAMICS 

To illustrate the influence of the initial conditions and 
spatial location of the atoms on the time development of 
the entanglement, we will use the solutions (|29|) and (|3Tj) 
to calculate the concurrence C(t) for some distances r 12 
and time t . For a fixed r 12 , we examine the time evolution 
of C(t) for the case where the atoms are first prepared in 
some pure state that can be specified by the Bloch vector 
components uq,vq and wq at time t — 0. In particular, 
we consider three different sets of initial pure states of 
the system corresponding to single-quantum excitations. 
In the first, the atom 1 is assumed to reside in its lower 
level |<7i) and the atom 2 in its upper level |e2), i.e. the 
initial conditions for the Bloch vector components are 
wq = 1, uq = vq = 0. In the second, the system is as- 
sumed prepared in a pure single-quantum superposition 
(entangled) state in which the two atoms oscillate in the 
opposite phase 

|*) = ^(|Si>|e 2 )+i|e 1 )| S2 )). (32) 

In this case, vq = 1, Uo = Wq = 0. In the third, the 
system is assumed prepared in a pure superposition state 
with the atoms oscillating with the same phase 

|^) = -L(| fll )|e 2 ) + | ei )| 52 )). (33) 

For this state, the initial conditions for the Bloch vector 
components are Uq = 1, Vq = u>o = 0. 

A. Preparation with too = 1> Uo = Vo = 

This limit corresponds to an initial condition in which 
the atom 1 is in the lower level and the atom 2 is in the 
upper level at time t = 0, i.e. in terms of the atomic den- 
sity matrix elements /922(0) = 1 and ^33(0) = ^23(0) = 
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P32(0) = 0. In this limit, the concurrence C(t) 
t = 0, and at later times is given by 



C(i) = e-^/l 



1 _ Ml sin 2 



1 



-at 



at 



(34) 



It is seen that the qualitative features of the transient 
development of the entanglement depend on whether or 
not 5i2 — 0. The concurrence sinusoidally varies with 
frequency a/2 whose the origin is in the dipole-dipole 
interaction and the presence of the detuning <5i2 ■ For the 
particular case of equal coupling constants when S12 = 
0, this result is reduced to that obtained previously in 
Refs. [IH[l6j]- The most obvious effect of having S12, i.e. 
unequal shifts of the atomic resonances, is that the the 
amplitude of the oscillating term is always less than 2. 
As a consequence, the entanglement may appear on the 
time scale twice as long as for the case of equal coupling 
constants. It is easy to see from Eq. (j34l) that for S12 = 0, 
the concurrence vanishes periodically at times 



io = nir/a, 



n= 0,1,2,. 



(35) 



whereas for 812 ^ 0, the concurrence behaves quite dif- 
ferently such that it vanishes only at times 



t$ = 2tg = 2nir/a, 



0,1,2, 



(36) 



i.e. the entanglement exists on the the time scale twice 
as long as for the case of equal coupling constants. 




FIG. 1: Transient evolution of the concurrence C(t) for 7 = 0, 
atom 1 located exactly at an antinode of the standing wave, 
and various locations 8r a of the atom 2 relative to an antin- 
ode of the standing wave: 8r a = (solid line), 8r a = 0.18A 
(dashed line). The system is initially in the state |2). The 
imperfect location Sr a = 0.18A, leads to the detuning 812 = 

2fii 2 . 

The features described above are easily seen in Fig. [TJ 
where we display the time evolution of C(t) for atom 1 
located exactly at an antinode of the standing wave and 
two different locations 8r a of the atom 2 relative to an 
antinode of the standing wave. According to Eq. (fPTj) . 
an imperfect or " nonideal" location of the atom 2 inside 
the cavity mode leads to a nonzero detuning 612- We 



assume a vanishing damping 7 = 0. This captures the 
essential dynamics of the system but makes no account- 
ing of the dissipative process during the transient evo- 
lution. We see oscillations in the transient evolution of 
the entanglement that follow the Rabi flopping of popu- 
lation back and forth between the atoms. In other words, 
this oscillation reflects nutation of the atomic populations 
which, in turn, can be associated with the precession of 
the Bloch vector B about the driving field vector f2 b with 
frequency a. 

The most interesting feature of the transient entangle- 
ment seen in Fig. [1] is that in the case of unequal cou- 
pling strengths, the initially unentangled system evolves 
into an entangled state, and remains in this state on the 
time scale twice as long as for the case of equal cou- 
pling strengths. This rather surprising result can be un- 
derstood in terms of spatial localization of the energy 
induced in the field by the initially excited atom. For 
equal coupling strengths the energy levels of the atoms 
are equally shifted by the amount <5i = 82 due to the in- 
teraction with the cavity mode. In this case the induced 
energy by the first atom oscillates with frequency a such 
that at the particular times t n — nix j a (n = 0, 1, 2, . . .) 
is fully absorbed by the second atom. Since at these 
times the energy is well localized in space as being com- 
pletely absorbed by the localized atoms, the entangle- 
ment, which results from an unlocalized energy, is zero. 
The situation changes when g\ ^ g 2 - According to (|2"7) . 
in this case the energy levels of the atoms are unequally 
shifted by the interaction with the cavity mode. Due 
to the frequency mismatch, the energy induced by the 
atom 1 is not fully absorbed by the atom 2, leading to 
a partial spatial derealization of the photon at discrete 
times t = nir/a, where n = 1,3,5,.... Consequently, at 
these times a partial entanglement is observed. The en- 
tanglement vanishes every time the excitation returns to 
its initial state, i.e. when it returns to atom 2. 



B. Preparation with vo = 1, uq 



■ Wo 







In this limit the initial population distributes equally be- 
tween the levels |2) and |3), i.e. in terms of the bare states 
density matrix elements are ^22(0) = P33(0) = P23(0) = 
P32(0) = 1/2. It then follows from Eqs. ^ and !(2TJ) 
that the concurrence C(t) = 1 at t = 0, and its time 
evolution is 



C(t) 



-it 



1 - 



4tt?2 . 2 



sin at. 



(37) 



We note immediately that the amplitude of the oscil- 
lating term that is equal to the population inversion is 
always less than unity provided S12 0. Thus, complete 
transfer of the population between the atoms cannot be 
achieved. In terms of the atomic excitation, none of the 
atoms can be completely inverted. This means that a 
part of the initial energy is always delocalized. Conse- 
quently, the system initially prepared in the entangled 
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state of the form (|32[) will remain entangled for all times. 
We emphasize that this feature arises from the presence 
of a non-zero detuning 5%2- That is the reason why the 
continuous in time entanglement is observed. Thus, this 
system may produce continuous in time atom-atom en- 
tanglement through having different values of the cou- 
pling constants g\ and g^- 




at 



FIG. 2: Concurrence C(t) as a function of time for 7 = 0, 
atom 1 located exactly at an antinode of the standing wave, 
and different 8r a : 5r a ~ (solid line), Sr a = 0.18A (dashed 
line). The system is initially in the state vo = 1, uo = U)q = 0. 

This behavior is illustrated in Fig. [2J where C(t) is plot- 
ted for various values of 812 and the initial superposition 
state ([521) . For t = the system is maximally entan- 
gled due to our choice of the initial state. Immediately 
afterwards, the entanglement begins to decrease because 
of the coherent oscillation of the atomic excitation. For 
the case S12 = 0, the system becomes unentangled peri- 
odically at the times t = nir/2, whereas for S12 ^ 0, the 
periodic minima (nodes) of C(t) are reduced in magni- 
tude resulting in a nonzero entanglement present for all 
times. 

The above analysis show that, rather surprisingly, im- 
perfect coupling of the atoms to the cavity mode may 
actually help generate continuous in time atom-atom en- 
tanglement, through unequal shifting of the atomic reso- 
nance frequencies. 



C. Preparation with uq = 1, vo = wo = 

If the system is initially prepared in the superposition 
state of the form (f33|) . the concurrence C(t) = 1 at t = 0, 
and the time evolution of the concurrence found from 
Eqs. ([2S]) and dSJ) is of the form 



C W =e^yi-^M^ S in 4 Q Q i). (38) 

It follows from Eq. ([55)) that in the absence of 612, i.e. 
when the atoms are in equivalent positions inside the 
cavity mode, the entanglement oscillation is completely 



suppressed. When 612 =/= the concurrence varies period- 
ically in time. The amplitude of the oscillating term that 
is equal to the population inversion is less than unity, un- 
less S12 — 2f2i2 and then the entanglement is completely 
quenched at the times t — mr/a. 



1 -r 
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FIG. 3: Transient evolution of the concurrence C(t) for 7 = 0, 
atom 1 located exactly at an antinode of the standing wave, 
and different 8r a - Sr a ~ (solid line), Sr a = 0.17A (dashed 
line). The system is initially in the state Uo = 1, vo = wo = 0. 

Figure [5J shows the influence of the detuning S%2 on the 
time evolution of C(t). It is evident that S12 significantly 
modifies the time evolution of C(t). For 612 = and 
vanishing damping, the entanglement is constant in time. 
Otherwise, when #12 7^ 0, the entanglement oscillates 
periodically and achieves the optimum value, C(t) = 1 
only at the particular times t — 2nn/a (n = 0, 1, 2, . . .). 
It can also vanish at the times t = nir/a. However, this 
happens only for a particular separation ryi for which 
<5i2 = 2f2i2. Thus, for positions of the atoms for which 
<5i2 7^ 2f2i2, entanglement is seen to occur over all times. 
This behavior of the entanglement is linked to the fact 
that with the pure entangled state preparation of the 
form ([3"3"| the Bloch vector and the dipole-dipole field 
vector are initially parallel. For 8\2 = 0, the Bloch vector 
B is effectively locked to the field vector, i.e. Bxfls = 0, 
and does not precess as it evolves on a time scale given 
by the spontaneous emission rate. In analogy to the spin 
locking effect, we may call this phenomenon as entangle- 
ment locking. When 5±2 7^ 0, the concurrence oscillates in 
time at frequencies a and a/2. In the Bloch picture this 
corresponds to the fact that for t > the Bloch vector B 
is no longer aligned along the VLb vector— precession of 
the Bloch vector translates into an oscillatory entangle- 
ment. 

Finally, let us examine in greater details the case 5\2 — 
2f2i2. We have seen that under this specific condition 
the concurrence exhibits interesting features. For exam- 
ple, in the case A, this is the value of 5\2 at which the 
entanglement achieves the optimum value with imper- 
fect matching. In the case B, the entanglement is always 
greater than 50%, and in the case C, the population is 
inverted for all times, and consequently the entanglement 
can be completely quenched at some discrete times. The 
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reason is that at this value of <5i2, the field vector Qb is 
in the uw plane, angle 9 = — 7r/4 from the north pole 
(w = 1) of the Bloch sphere. The Bloch vector processes 
about a cone whose opening angle 9 depends on the ini- 
tial conditions. In the case A, the the Bloch vector makes 
a constant angle 9 = ir/4 with Hb, whereas in the case C, 
it makes a constant angle 9 = 7r/4 with — fls- Thus, in 
these two cases, the Bloch vector rotates in a quarter of 
the Bloch sphere such that it can reach one of the poles, 
w = ±1. It therefore appears that in the case A, the 
Bloch vector when processing about f2_e it varies from 
w = to w — 1, i.e. regularly reaches the north pole, 
but in the case B, it varies from w — to w = — 1, i.e. 
regularly reaches the south pole. The case C is different, 
the Bloch vector makes an angle 9 — tt/2 with CIb, so 
that it processes about Qb in such a way that it rotates 
from w = — 1/ V2 to w = 1/ \[2 such that it never achieve 
the poles w = ±1. 



V. ENTANGLEMENT DIFFRACTION 

One of the more interesting aspects of the transient en- 
tanglement demonstrated in the previous section is its 
dependence on the difference of the Stark shifts of the 
atomic transition frequencies. This difference arises from 
an imperfect coupling of the atoms to the cavity mode 
that is a consequence of noncquivalcnt positions of the 
atoms inside the cavity mode. We proceed here to present 
more detailed studies of the sensitivity of the transient 
entanglement to position of the atoms inside a standing- 
wave cavity mode. We choose the initial conditions to 
be wo = l,ito = vo = 0. When Eqs. © and (fTTj) are 
used in Eq. (j3~4"|) , we readily find that the variation of the 
concurrence with position of the atoms in the standing 
wave is given by 



C(ria) 



,-rv 



'sine? 



+ 



sin i<i 



t 4 sin 4 fcr 12 



| cosfcri 2 |, (39) 



where d = (1 + cos 2 kri2)r/2 and we have introduced a 
scaled time variable 



A 

and the dimensionless damping rate 



(40) 



(41) 



both measured in units of 2<? 2 /A which is always assumed 
to be much smaller than unity. 

The concurrence (|39[) exhibits an interesting modulation 
of the amplitude of the harmonic oscillation. One could 



naively think that a variation of the concurrence with 
Tyi should reveal the cosine form of the cavity mode 
function. However, Eq. (|39|) shows that the concurrence 
is not a simple cosine function of r\2- It is given by 
the product of two terms, one the absolute value of the 
cavity mode function | cos kr\2 | and the other the time- 
and position-dependent diffraction structure. That is, 
the amplitude of the standing wave cavity mode is in 
the form of position and time dependent diffraction pat- 
tern. Only at very early times (r <C 1), the oscillations 
are not modulated by the diffraction pattern and the 
concurrence reduces to |cos/cri2|, but for longer times, 
C(r 12) may vary slower or faster than the cosine func- 
tions. Within the diffraction structure itself, the magni- 
tude of the concurrence exhibits a succession of modes 
and of antinodes. As a consequence, the entanglement 
may be completely quenched even for locations of the 
atom 2 close to an antinode of the cavity mode, and alter- 
natively may achieve its optimum value even for locations 
of the atom 2 close to a node of the cavity mode. It is 
easy to show that C(r 12) vanishes periodically whenever 
cos [(l + cos 2 fcri 2 ) t/2] = 1, i.e. for discrete times 



Awk 



1 + cos 2 kr\2 



n = 1,2,3, . 



(42) 



We may establish a relation between the number of ze- 
ros in C(r 12) and the time r by an elementary argument. 
Since cos kryi varies between zero and one, the shortest 
time at which C(r 12 ) achieves at least one zero is r = 27r. 
Thus, there are no zeros in C(tyi) for locations of the 
atom between the successive nodes of the cavity mode 
function, if r < 2tt. For a given r the number of zeros 
is limited by the fact that cos 2 kr 12 < 1. This imposes 
lower and upper limit on n in Eq. (|42|) . The upper limit 
is important, since it determines the width of the region 
about antinodes of the cavity function where the opti- 
mum entanglement occurs. It follows from Eq. (|42p that 
the criterion for vanishing entanglement is satisfied for 



T 

n < 

~ 2tt 



(43) 



The largest n satisfying this inequality corresponds to 
the largest value of cos 2 kr\2 , and therefore determines 
a node that is the closest to the antinode of the cavity 
mode function, so that it determines the width of the 
main peak of the diffraction pattern. 
More interesting is a possibility of obtaining the optimum 
entanglement when the atom 2 is located between a node 
and a successive antinode of the cavity mode. The op- 
timum entanglement occurs at the locations of the most 
intense maxima of the concurrence. However, the loca- 
tions of the maxima are not given by a simple relation. 
For a given t, the concurrence (|39[) achieves the optimum 
value C(r 12) = 1 whenever 



cos 



- (1 + cos 2 kr 12 ) t 



sin kr\2 
2 cos kr\2 



(44) 
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This is not a simple relation, and we solve this equation 
graphically as follows. Introducing the notation 



q(n 2 ) 



i 



(l + cos 2 krw) t 



sin 2 fcri2 
2 cos kr\2 



(45) 



we find solutions of the equation p(ri 2 ) = q(r 12) by plot- 
ting separately p{ri%) and q{r\2). The functions p(rvi) 
and q{r\2) are shown in Fig. 3J The intersection points 








0.4 „ 0.6 

r 12 /X 



FIG. 5: Concurrence C{r\_2) as a function of the separation 
between the atoms when 7 = 0, the system is initially in 
the lower level, wo = 1, ito = vo = 0, and (a) r = it/2, (b) 
t = 97r/2, (c) r = 27tt/2. 



entanglement oscillates with position faster than the co- 
sine function, and the oscillations are more dramatic for 
larger times. 



FIG. 4: The parameters q(r\2) (solid line) and p(r\2) plotted 
as a function of ri2 for different times r: ri = tv/2 (dashed 
line), (b) r 5 = 9ir/2 (dashed-dotted line), (c) no = 277r/2 
(dotted line). The choice of the particular values of r corre- 
sponds to optimum entanglement observed for the idealized 
case of gi = gi. 

of the two curves give the solutions of the equation (|4"4")) . 
At these points the system attains the optimum entan- 
glement. We see from the figure that the equation (144]) 
is satisfied only for discrete values of tyi. The number 
of solutions, which gives us the number of the optimum 
that the concurrence may achieve, depends on time r. 
Rather than examine the situation at all times we will 
look only at the particular times r„ = mr/2, correspond- 
ing to the evolution intervals the optimal entanglement 
is obtained for the idealized case of g\ = g2- For r = tt/2 
there is only one solution corresponding to the position of 
the atom 2 precisely at the antinode of the cavity mode 
function. For a longer time r = 97r/2 there are three 
solutions and the number of solutions increases with r. 
Figure [5] shows C(t n ) as a function of r^/A for differ- 
ent times r. For a short time the entanglement is seen 
to occur over a wide range of positions centered about 
the antinodes of the cavity mode. The concurrence is a 
bell-shaped function of position without any oscillation. 
As time progresses, oscillations appear and consequently 
the region of r% 2 where the optimum entanglement oc- 
curs, becomes narrower. The evolution of C(r 12) tends 
to become increasingly oscillatory with r±2 as time in- 
creases, and the optimum entanglement occurs in a still 
more restricted range of r%2. As a result, the atom-atom 




FIG. 6: Concurrence C(ri2) as a function of the separation 
between the atoms when 7 = 0, the system is initially in the 
lower level, too = 1, uo = Vo = 0, and (a) r = n, (b) r = 10tv, 
(c) r = 28tt. 

In Fig. [6l we plot C{t n ) as a function of r^/A for three 
values of time r corresponding to the evolution intervals 
at which the entanglement is quenched for the idealized 
case of gi = 92- Here we observe, that the concurrence 
actually evolves with the position of the atom leading to 
the appearance of what we may call an inverse diffraction 
pattern. Note that the concurrence vanishes for locations 
of the atom precisely at the antinodes of the cavity mode, 
and may achieve its optimum value at locations of the 
atom close to the nodes of the cavity mode. 
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This effect can be understood as a consequence of the 
uncertainty relation between the evolution time and en- 
ergy [23|, [13] . For the increasing time the uncertainty of 
the energy decreases which means that the energy be- 
comes more localized. The increase in the localization of 
the energy results in a degradation of the entanglement. 
In other words, with increasing time, one can in princi- 
ple obtain more information about the localization of the 
atoms inside the cavity mode. 

VI. SUMMARY 

We have investigated the process involved in the entan- 
glement creation between two distant atoms coupled to 
a single- mode cavity field. Unlike previous publications, 
we have included a possible variation of the coupling con- 
stant g(r) with the location of the atoms in a standing- 
wave cavity mode. We have found that the entanglement 
creation in a complex two-qubit system can be modeled 
in terms of the coherent dynamics of a simple single- 
qubit system driven by a coherent field. Effectively, we 
have shown how to engineer two coupled qubits whose 
the dynamics are analogous to that of a driven single two- 
level system. We have obtained analytical expressions for 
the concurrence and have shown some new properties of 
the entanglement that are not met in the idealized case 
of equal coupling constants appear for unequal or "im- 
perfect" coupling constants. In particular, the degree of 
entanglement and the time it takes for the concurrence 
to reach its optimum value is a sensitive function of the 
position of the atoms inside the cavity mode. Charac- 



terizing the system by the Bloch vector components, we 
have examined the parameter ranges in which entangle- 
ment can take place for all times. We have demonstrated 
that a spatial variation of the coupling constant affects 
localization of the energy induced in the field by interact- 
ing atoms that leads to a long-lived entanglement. The 
consequence of this imperfection is that under certain ini- 
tial conditions, an initially entangled system may remain 
entanglement for all times. 

Finally, we have shown that the variation of the concur- 
rence with the position of the atoms is that of the cavity 
mode function multiplied by a time-dependent diffrac- 
tion pattern. The diffraction formula shows explicitly 
the trend of the modification of the entanglement with 
the localization of the atoms when the observation time 
increases. For a short time the entanglement is seen 
to occur over a wide range of positions centered about 
the antinodes of the cavity mode. As time progresses, 
oscillations appear and consequently the spatial region 
where the optimum entanglement occurs, becomes 
narrower. This effect has been explained in terms of 
the quantum property of complementarity, which is 
manifested as a tradeoff between the knowledge of 
energy of the exchanged photon versus the evolution 
time of the system. 
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